Van der Waals forces and spatial dispersion 
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A version of the Green's functions theory of the Van der Waals forces which can be conveniently 
used in the presence of spatial dispersion is presented. The theory is based on the fluctuation- 
dissipation theorem and is valid for interacting bodies, separated by vacuum. Objections against 
theories acounting for the spatial dispersion are discussed. 
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Introduction. It was recently discovered that effects 
of the spatial dispersion are quite important for calcula- 
tions of the Van der Waals forces [ll between conducting 
bodies at finite temperatures (see [2, 0, 0, 0, 0]). On the 
other hand, an opinion has been expressed in several pa- 
pers that the theory of the forces in the present form is 
not valid in the presence of this dispersion (see 0,0,0]). 
The goal of this paper is to prove rigorously that the the- 
ory of the Van der Waals forces, as it was formulated in 
the terms of Green's functions in papers 13, II | (DLP 
below, see also [13), can be used for calculation of the 
forces for bodies separated by vacuum, even if the spatial 
dispersion is important for electromagnetic properties of 
the bodies. 

To my understanding, the doubts about the validity 
of the theory in the presence of the spatial dispersion is 
based on a remark in a paper by Barash and Ginzburg 
[ljj . In this paper, the contribution to the free energy of 
the long- wave fluctuations of the electromagnetic field is 
presented in the form 
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where V n is the Matsubara Green's function of the sys- 
tem under considerations and T> n (£) and II„ (£) are cor- 
respondingly the Green's function and the polarization 
operator of an auxiliary system in which interaction with 
the long- wave field is defined by the " charge" VCe instead 
the actual charge e (see (l3| . Eq. (54)). The prime sign 
means that the n = term is taken with the coefficient 
1/2. The first term in JT]) corresponds to the DLP theory, 
while the second gives an additional contribution. The 
authors noticed that in the presence of significant spa- 
tial dispersion the second term cannot be neglected and 
consequently the free energy cannot be expressed only in 
terms of the Green's functions of the actual system. 

However, in my opinion this statement is related to 
the general case, when bodies are embedded in a liq- 
uid with spatial dispersion. Results of the DLP theory 
are indeed difficult to generalize to this case due to the 



non-local nature of the stress tensor in a such liquid. In 
the important case when the interacting bodies are sep- 
arated by vacuum this difficulty does not arise. Indeed, 
forces acting between any bodies separated by vacuum 
can be calculated by averaging of the vacuum Maxwell 
stress tensor. Corresponding quadratic combinations of 
the field strengths can in turn be expressed in terms of 
the retarded Green's function of the field using the exact 
fluctuation-dissipation theorem (see [12]). This Green's 
function, of course, must be calculated taking into ac- 
count the spatial dispersion in the bodies, when it is 
important. This programme wil be realized below. It 
is convenient to use an approach which is different both 
from nj, 11 1 and [T3|. (The basic idea of the derivation 
was expressed in short in Q.) 

Maxwell stress tensor. The force acting in vacuum 
on a body due to electromagnetic fluctuations is equal 
to 14] 
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where the integration is performed over any surface en- 
closing only the given body and o~i k is the average value of 
the Maxwell stress tensor of electromagnetic fluctuations 
in vacuum: 

(E 2 ) - (E t E k ) (B 2 ) r (BtB h ) , s 



8tt 



4tt 



8tt 



4tt 



Here we introduced notations for the averages from op- 
erator expressions: 



(EiE k ) = {Ei(r,t)E k (r,t) + E k (r,t)Ei(r,t)j /2 ,(4) 

(E 2 ) = {EiEi) 

and 

{B l B k ) = (Bi(r,t)B k (r,t) + B k (r,t)Bi(r,t)} /2 ,(5) 

(B 2 ) = (BiBt) . 

I consider only the equilibrium situation and averaging 
is taken using the Gibbs statistics. Thus the problem 
reduces to calculation of the average of quadratic com- 
binations of field strength operators. Notice, that this 
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method of calculation of forces can be applied only in 
vacuum, because there is in general no equation for the 
stress tensor of non-static fields in an arbitrary medium. 

Fluctuation-dissipation theorem. To define these 
averages one can use the exact fluctuation-dissipation 
theorem (FDT). To apply this theorem one must con- 
sider correlation functions of the fields in different points 
in space and different moments of time: 

^(r,r';r) (6) 
E { (r, t + r)E k (r' , t) + E k (r', t)E t (r , t + r) \ /2 , 



where we took into account the uniformity with respect 
to time. The magnetic function S 3 , (r, r'; r) is defined in 
an analogous way. We need also the correlation function 
of the components of the 4-vector potential. Following 
fioj ] we will use the gauge with the scalar potential (p = 0. 
Then 



f OA ~ - 
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and Sf k , SR can be expressed as (we put c = h = 1 
below) 



Sg (r, r'; r) = (curl) is (curl') fc; (r, r'; r) 



(8) 



in the terms of the correlation function of the components 
of A: 



S-t(r,r';r) (9) 
ii(r, t + r)i fe (r', t) + i fe (r', t)li(r, t + r)\ /2 . 



Expand now the function S^. in the Fourier integral with 
respect to time: 
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The central point of the proof is that the fluctuation- 
dissipation theorem allows us to express S-j, (r, r';ui) in 
the terms of the retarded Green's function of the electro- 
magnetic field. Indeed, the interaction of the field with 
the current is given by the equation 
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Let us consider the vector-potential which is induced by 
a classical external current with the density j ext (r, t) — 
[jS"* ( r ) e~ luJt + (r) e iu)t ] /2. In the linear approxima- 
tion the induced potential can be expressed in the terms 
of the retarded Green's function of the vector-potential: 



(12) 



where, according the Kubo equation, 



(13) 
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Then according to the fluctuation-dissipation theorem 
(see [I*!, sections 122-125 and 12|, sections 75-76), 
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If the bodies are not magnetoactive, the £>-function 
satisfies the Onsager symmetry relation D R k (r, r';u>) — 
(r', r; w) ■ Then we find finally 

S& (r, r» = - coth (-|^) ImDffe (r, r'; W ) . (15) 

This equation together with ([8]) solves the problem of 
calculating electromagnetic fluctuations in the terms of 
D R , One has 
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Equation (| 13[) is not very useful for calculating the 
function Dg. However, D R can be calculated directly 
using equation (|12p . Indeed, let us place a unit point- like 
current directed in the /-direction at the point r' = ro, i. 
e. put 



x [(curl) j;s (curl') fcI D*(r,r i; a,)] 
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(18) 



Using a proper theory of electrodynamic properties of the 
interacting bodies, one must calculate the linear response 
A^j (r, r ) for the source Then according to (TR 

wc get 



Dg (r,r ;o;) = ^ (r,r ) 



(19) 



It is important that it is enough for the calculation of the 
tensor to consider the source placed in vacuum between 
the bodies. This, for example, allows in a simple plain 
geometry to express in the terms of the reflection 
amplitude on the surfaces of the bodies. Expression (TIT))) 
for Dg diverges at r — > tq. This divergency can be elim- 
inated by subtracting the response A™£ (r, r ) for the 
source (fT5)). calculated for free space. 
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Notice that according to (fl6|) - (flT|) the stress tensor 
can be presented in the form 

^=^f^^)af M ^ (20) 

where <tt| (u>) is an analytical function having no singular- 
ities in the upper half-plane u>. Displacing the contour of 
integration to the imaginary axes of u, as in the original 
Lifshitz paper [12], and taking the residues in the poles of 
the coth, we reduce ([20]) to the Matsubara form: 

Conclusion. The general theory which I present here 
does not give, of course, a magic prescription for solution 
of all problems of Van der Waals forces in the presence 
of spatial dispersion. Calculations of the D R function 
as a response to a point source can be difficult and in 
any case demands a sort of microscopic theory. How- 
ever, I believe that reducing the problem to the exact 
fluctuation-dissipation theorem gives a solid foundation 
for different theoretical approaches. 

It must be noted however that in a recent paper Q the 
authors express doubts about the validity of the Lifshitz 
theory based on this theorem for materials with finite 
d.c. conductivity. Because FDT is a rigorous theorem of 
quantum statistical mechanics, such doubts must have a 
very serious foundation. On the contrary, the considera- 
tion of the authors are, in my opinion, quite superficial. 
One can read in [§] (I omit references): 

"Lifshitz derived his famous formulae under the con- 
dition of thermal equilibrium. This means that not only 
T = const, but also all irreversible processes connected 
with the dissipation of energy into heat have already been 
terminated... The Drude- like dielectric function ... is de- 
rived from the Maxwell equations with a real drift current 
of conduction electrons j = ooE initiated by the exter- 
nal electric field E ... The drift current is an irreversible 
process which brings a system out of thermal equilib- 
rium. ... The real current leads to Joule's heating of the 
Casimir plates (Ohmic losses) ... To preserve the tem- 
perature constant, one should admit that there exists an 
unidirectional flux of heat from the medium to the heat 
reservoir ... Such interactions between a system and a 
heat reservoir are prohibited by the definition of thermal 
equilibrium. Although the screening and diffusion effects 
really occur in an external electric field, they are also re- 
lated to physical situations out of thermal equilibrium. 
The reason is that the diffusion current is determined by 
a nonzero gradient of charge carrier density, whereas for 
homogeneous systems in thermal equilibrium the charge 
carrier density must be homogeneous." 

I believe that these considerations are wrong. Of 
course, one can say that the fluctuating electric field 
heats a body. However, in equilibrium this heating is 



compensated by emission of radiation by the body. This 
exact compensation is ensured by the detailed balance 
principle. The sentence about termination of dissipation 
of energy is particularly odd in relation to the fluctuation- 
dissipation theorem, which just connects the energy dis- 
sipation with fluctuations. It is not clear, by the way, 
why the authors worry only about the Ohmic losses. At 
finite frequencies all real materials dissipate energy. 

The statement that the screening and diffusion effects 
are related to physical situations out of thermal equilib- 
rium is also wrong. It is well-known that the Boltzmann 
distribution in the electric field, which was used in Q 
for describing screening, is an equilibrium distribution. 
Actually, in equilibrium the diffusion current is compen- 
sated by the mobility of carriers due to an electric field. 

The authors of [9( claim also that "for dielectrics whose 
charge carrier density is temperature-independent (for 
such materials conductivity goes to zero with T not due 
to the vanishing n but due to the vanishing mobility of 
the charge carriers) the generalization of the Lifshitz the- 
ory taking into account the screening effects is shown 
to violate the Nernst theorem" and is consequently "in 
contradiction with thermodynamics" . I believe that this 
statement is a result of a pure misunderstanding. The 
materials under discussion are amorphous glass-like dis- 
ordered bodies. Conductivity goes to zero with T in such 
materials due to localization of the charge carriers just 
because of the disorder. The point is that the Nernst 
theorem is not valid for these disordered bodies. It is 
well known that they have a big finite entropy at zero 
temperature. Localized carriers also contribute to this 
residual entropy and the calculation of a small correction 
to its value due to the Van der Waals interaction scarcely 
has a physical meaning. Of course, the existence of disor- 
dered bodies at T = itself does not contradict statistical 
mechanics. They are simply not at an equilibrium state 
at low temperatures due to a very long relaxation time. 
However, the role for our phenomena of the long relax- 
ation time in these bodies is not clear at present and is 
worth of careful investigation. 

I thank Yu. Barash, V. Svetovoy and E. Taylor for 
useful discussions. 
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